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TECHNICAL NOTE 3864- 

THEORETICAL CALCULATION OF THE POWER SPECTRA 
OF THE ROLLING AND YAWING MOMENTS ON A 
WING IN RANDOM TURBULENCE 
By John M. Eggleston and Franklin W. Diederich 

SUMMARY 


The correlation functions and power spectra of the rolling and yawing 
moments on an airplane wing due to the three components of continuous 
random turbulence are calculated. The rolling moments due to the longi- 
tudinal (horizontal) and normal (vertical) components depend on the span- 
wise distributions of instantaneous gust intensity, which are taken into 
account by using the inherent properties of symmetry of isotropic turbu- 
lence. The results consist of expressions for correlation functions or 
spectra of the rolling moment in terms of the point correlation functions 
of the two components of turbulence. 

Specific numerical calculations are made for a pair of correlation 
functions given by simple analytic expressions which fit available experi- 
mental data quite well. Calculations are made for four lift distributions. 
Comparison is made with the results of previous analyses which assumed 
random turbulence along the flight path and linear variations of gust 
velocity across the span. 

The rolling moment due to lateral (side) gusts, which is small, is 
expressed in terms of the instantaneous value of the gust near the center 
line of the fuselage, so that the effect of spanwise variation in gust 
Intensity is ignored. The yawing moments are considered to be propor- 
tional to the rolling moments with the constants of proportionality given 
by simple aerodynamic relations. 


INTRODUCTION 


The gust velocities acting on an airplane flying through turbulent 
air are functions of position or time known only in a statistical sense. 
Consequently, aerodynamic forces and moments produced by the lifting 
surfaces of the airplane can be known only in a statistical sense. If 
the statistical characteristics of the turbulence are assumed to be 
invariant with position along the flight path, flight through turbulent 
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air may be considered to be a stationary random process and the mathe- 
matical techniques developed for such processes (see ref. 1, for instance) 
may then be used in this problem. 

This approach has been adopted in many papers on this subject, among 
them references 2 and 3* Inasmuch as in these papers the motions and 
forces associated with the longitudinal degrees of freedom were of primary 
interest, the assumption was made, implicitly, that the gust intensity is 
uniform along the span at any instant. However, for the problem of ana- 
lyzing the motions and forces associated with the lateral degrees of free- 
dom, this assumption is inadequate, inasmuch as it implies that the verti- 
cal. and horizontal gusts produce zero rolling and yawing moments on the 
wing. This problem has been treated in references 4- and 5 and elsewhere 
on the basis of the assumption that at any instant the guBt intensity 
varies linearly across the span. 

A fundamental method of accounting for the lift on a wing due to 
random variations of the gust velocities in both the flight-path and the 
spanwise directions is given in reference 6 for the longitudinal response 
of bn airplane in atmospheric turbulence. The approach is based on the 
assumption that the turbulence is axisymmetric (according to ref. 7 ) > so 
that, at any arbitrary time or position in the turbulence, the statistical 
characteristics of the turbulence encountered by an airplane do not depend 
on the heading of the airplane. On the basis of this assumption, the 
variation of gust intensity across the span can be related to the varia- 
tion of the gust intensity along the flight path. 

In the present paper the approach of reference 6 Is extended to the 
calculation of the rolling and yawing moments on a wing due directly to 
vertical gusts, longitudinal gusts (hereinafter referred to as horizontal 
gusts), and lateral or side gusts. These moments are required as a first 
step in calculating the motions of a complete airplane in atmospheric 
turbulence; the moments due to the motions caused by these input moments 
can be calculated by conventional methodB and will not be considered 
herein. 

In the first part of the paper, a theoretical analysis is made 
defining the power spectra of the rolling and yawing moments of a wing 
in terms of the statistical characteristics of the atmospheric gust veloc- 
ities. By using. an analytical expression to define these characteristics, 
a numerical solution of the lateral moments is presented in the last 
part of the paper. 
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E(k), K(k) 


f 

F 

g 

G 

h 

h' 

I 

1 -^ 

k 


k* 

*0* K X 

<V A/ 

* 0 , K 1 


wing span 
wing chord 

wing mean aerodynamic chord 

complete elliptic integrals of the second and first kind, 
respectively, of modulus k 

longitudinal correlation function for isotropic turbulence 
Fourier transform of f 

lateral correlation function for isotropic turbulence 
Fourier transform of g 
indicial-response function of time only 
indie ial-response function of time and displacement 
Fourier transform of two-dimensional correlation function 

modulus of elliptic integrals, — ~ — 

2 + T) 

reduced frequency, oi/U 

modified Bessel functions of the second kind 
incomplete modified Bessel functions of the second kind 


l 

L 

Mx 

P 


section lift 

integral scale of turbulence 
rolling moment 
rolling velocity 
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q. 

r 


S 

t 

U 

Ur 

v 

UgJ Vg, Wg 


dynamic pressure 

yawing velocity (used only in stability derivatives); 
linear displacement between any two points 

wing area 

time 

mean forward velocity 

displacement along the flight path 

component of airplane velocity along positive Y-axis 

three components of gust velocity (see fig. 1(a)) 


X, Y, Z 
x 

y 


reference axes (see fig. 1(a)) 
chordwise distance 
spanwise distance 


Ar - y 2 “ Yi 


y* 

a 

P ' ~ b/L 
7 

r 

A = Ut/L 

t] 

P 

T 

03 


nondimensional spanwise coordinate, 

b/2 

angle of attack, radians 

span influence function 
integral weighting function 

dummy variable of integration, - y-j* 

atmospheric ■ density 
dummy variable of time 
circular frequency, 2it/Period 
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C2 


rolling-moment coefficient. 


Rolling moment 
qSb 


Cn 


yawing-moment coefficient. 


Yawing moment 
qSb 


c n p 




bCi 



2U 


Cn r 


6Cn 

2U 


i 





2U 





8C i 
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correlation function 
power spectral density 
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Subscripts : 

o trim, value 

g gust component 

A bar over a quantity denotes the mean value of the quantity. The 
absolute value of a quantity is denoted by | [ . 


THEORETICAL ANALYSIS 


Preliminary Considerations 

In this section expressions are derived for the power spectra of 
the rolling and yawing moments of an unswept airplane wing or thin lifting 
surface of arbitrary plan form due to flight through random atmospheric 
turbulence. Essentially, the procedure consists of expressing the rolling 
moment at any arbitrary position along the flight path in terms of the 
gust velocity at that position, establishing the correlation function 
between the rolling moments at any two points along the flight path, and 
transforming this correlation function into an expression for the power 
spectral density. The power spectrum of the yawing moment is then related 
to that of the rolling moment through simple aerodynamic relationships. 

Assumptions . - The following assumptions are made in the analysis: 

(1) The turbulence is homogeneous and isotropic; that is, the statis- 
tical characteristics of the turbulence are invariant under a translation 
or rotation of the space axes (although the results obtained for the verti- 
cal component of turbulence require only the somewhat less restricting 
assumption of axisymmetry) . 

(2) Time correlations are equivalent to space correlations along 

the flight path - an assumption usually referred to as Taylor's hypothesis. 
(See ref. 7 .) 

(3) The chordwise penetration factor (the indicial-response influence 
function) for the rolling and yawing moments can be expressed as a product 
of a function of distance along the flight path (or time) only and distance 
along the span only. 

(4) The wing considered herein is relatively rigid and, as a result 
of the turbulent velocities, performs small motions about a mean steady 
flight condition. 

The implication of these assumptions and the limitations they impose 
on the results of the analysis are discussed in a subsequent section of 
the paper. 
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Coordinate system and, gust components . - The system of axes and the 
local velocity field relative to the lifting surface are shown in fig- 
ure l(a) . The velocity at each point in the field is resolved into com- 
ponents lying in the three planes of an orthogonal set of axes, the X-axis 
of which is tangent at every point to the flight path. Throughout this 
paper these three components are designated as follows: The component 

alined with the X-axis is referred to as the horizontal gust Ug; the 

component alined with the Y-axis is referred to as the side gust v g ; 
and the component alined with the Z-axis is referred to as the vertical 

gUSt Wg. 

As the wing moves through the local velocity field, the random vari- 
ations in the horizontal and vertical gust components are defined both 
in the flight-path direction and in the spanwise direction at every posi- 
tion along the flight path. Random variations of these gust components 
across the chord are taken into account by indicial-response functions 
and, hence, need not be considered separately. 

The side gust component of the gust velocity field is treated in 
only a limited manner. Neither the chordwise nor the spanwise variations 
of Vg are considered along the flight path; rather, v g is assumed to 

act on the wing as a point velocity with a variation only along the 
flight-path direction. Contemporary aircraft exhibit such wide varia- 
tions in distribution of dihedral across the span that it is doubtful 
that a generalized analysis could be utilized. The point or centroid 
analysis should be fairly accurate when the dihedral distribution is 
predominant over only a small section of the span near the fuselage. 

Such a distribution is exhibited by an unswept wing with zero geometric 
dihedral mounted very high or low on a fuselage. For a wing with zero 
aerodynamic dihedral, this component could be neglected completely. 

Definition of gust correlation functions .- In order to define random 
variations of the gust velocities both along the flight path and across 
the span of the wing as it moves through the turbulence, it is necessary 
to define the correlation between any two velocities in the gust field 
through which the wing passes. The space correlation function of a veloc- 
ity u is defined in terms of the distance r as 



( 1 ) 


Von Karman and Howarth (ref. 8) have shown that, in homogeneous isotropic 
turbulence, the correlation between two velocity vectors a distance r 



8 


NACA TN 3864 


apart can be defined in terms of two scalar functions f(r) and g(r) 
and that this relationship is invariant with respect to rotation and 
reflection of the coordinate axes., These one -dimensional correlation 
functions relate the paired velocity components obtained by resolving 
the velocity vector at any two points a distance r apart into two parts 
The pair lying along the straight-line path between the points are known 
as the longitudinal components and the pair normal to the straight-line 
path are known as the lateral components . These two pairs of components 
are pictorially shown in figure l(b). Such velocity components may be 
measured in wind tunnels downstream of a grid mesh. (See ref. 9-) 

In reference 8, it is further shown that these correlation functions 
are interrelated by the differential equation 


1 Mil + f(r ) , g (r) (2) 

2 dr 

By defining the variable r in the coordinate system of this paper and 
using the correlation tensor of reference 8, a two-dimensional analysis 
of the turbulence as it affects the wing may be made in terms of f(r) 
and g(r) . The variable in the correlation functions of the horizontal 
and vertical gust components in the two-dimensional XY -plane of the wing 
is given simply by 

r = )/(£x) 2 + (Ay) 2 = ^(Ut) 2 + (Ay) 2 (3) • 

The correlation function of the horizontal gust components , as 
derived from the correlation tensor of reference 8, Is defined in terms 
of x- and y-components of the present analysis by the formula 


$ Ug (^Asr) - ug <j 


r 


(^) a f 

\j(Ax ) 2 + (Ay) 2 

(£x ) 2 + (Ay ) 2 

Jf 


(Ay) ! 


(/he) 2 + (ay) 2 


s 


^j(£x ) 2 + (dy)‘ 


(4) 


The relationship between the components is shown schematic ally in fig- 
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In a like manner, the correlation function of the vertical gust 
components affecting the wing, given in terms of the mean-square value 

of the vertical gust velocity Wg^, may be seen to be simply 


$ Wg (£x,4y) 


2 

w g s 




Zsx) 2 + (Zy)‘ 


( 5 ) 


For the case of side gusts acting on a wing, the correlation func- 
tion would be defined in terms of Zac and Zy if the spanwise correla- 
tion were considered. (See fig. 2(b).) Inasmuch as the side gust is 
considered to act only at a point on the span, Zy is zero, and the 
correlation function for the side gust in terms of its mean-square value 
becomes 


$ Vg (Zcc) = v g 2 g(Zsc) (6) 

Although the mean-square value of each of the three gust components 
is given separate identity, under the assumption of isotropy 


Ug2 = v g 2 = w g 2 (7) 

With the gust-velocity correlation functions thus defined, the 
forces and moments due to antisymmetric components of the gust-velocity 
field acting on a wing passing through that field may be derived in terms 
of these correlation functions. 


Rolling Moment Due to Gusts 

Vertical gusts .- The instantaneous wing rolling moment due to vertical 
gusts can be written in terms of an indicial-response influence func- 
tion h 1 (t,y) as 


*b/2 


Mx(t) = J J ^'(t^y) Vgjt-t^y) dy dt x 


( 8 ) 


According to assumption (3) of the section entitled "Preliminary Con- 
siderations" (see also the argument presented in ref. 6), the func- 
tion h'(t,y) can be expressed in the form 
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h' (ti,y) = h ( t i) 7(y) (9) 

where y(y) is a steady-state span, influence function and h^tq) con- 
tains the unsteady- lift effects. The rolling moment can then also be 
written as 


%(t) = J h(t x ) dt x J* r(y) Wg^b-ti,yj dy (10) 


If the correlation function for the rolling moment is calculated 
from equation ( 10 ) and a power spectrum for the rolling moment is obtained 
by taking the Fourier transform of the correlation function, the resulting 
expression may be shown to consist of a product of two functions: One 

function is the result obtained from quasi-steady considerations alone, 
and the other is the absolute squared value of the unsteady- lift function 
for sinusoidal gust penetration such as that given by Sears in refer- 
ence 10. Consequently, consideration will be confined to an analysis 
using quasi-steady expressions for the rolling moment; that is, the lag 
in buildup of lift across the chord of the wing due to the gusts is not 
included. 

In quasi-steady flow, the rolling moment of a wing due to a variable 
angle-of -attack distribution across the span is given by 


M X = qSbCj 

■C 


where section lift 


( 11 ) 


l(y) = Cj(y) q c(y) 


and local angle of attack due to gusts 
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How, one theorem of linearized airfoil theory states that the lift (or 
rolling moment) on a wing due to an arbitrary spanwise angle-of- attack 
distribution is equal to the integral over the entire wing of the product 
of the spanwise lift distribution due to a unit constant (or linearly 
varying) angle of attack and the given arbitrary angle -of- attack distri- 
bution. Hence, the rolling moment is also given by 


qSb Ci 



r- -i a?=y 

Uy)J ag(y) ciy 


( 12 ) 


This theorem is valid not only in steady but also in indicial flow. 

(See reciprocity theorems of ref. 11.) 

When the indicated substitutions are made, the rolling-moment coef- 
ficient along the flight path is 


C}(x) 


•c 7 

ob/2 

Cq(y) c(y) 

b^ 

^-b/2 

-Clp 8 


cc=y 


w g (x,y) 


U 


<3y 


pl 


4u 


J 7(y*) w g (x,y*) dy* 


(15) 


where y* = -*4_ and the steady-state lift distribution 
b/2 


y(y* ) = 



(y*) 1 c (y*) 


ct=y* 


( 14 ) 


pertains to a linear antisymmetric angle of attack across the span. It 
may be seen that, by virtue of its definition, y(y*) must satisfy the 
relation 


7 (y*) y* dy* = 2 

0 


(15) 
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Horizontal gusts .- In analogy to the analysis of the preceding sec- 
tion, consideration will be confined to the quasi-steady case. When 
stability axes are. used, a change in forward velocity at any spanwise 
station increases the magnitude but does not change the direction of the 
lift and drag vectors. Thus, the horizontal-gust contribution to the 
dynamic pressure is 


Aq(y) = 



U 2 


- \ p( u s 2 + atog) 


— 2qQ 


ug(y) 

u 


under the assumption that Ug « U. When this linearized approximation 

is used, the lift on each section is proportional to the local angle of 
attack : 


Z(w g ) w g 

The rolling-moment coefficient due to horizontal-gust velocities is 
thereby defined as 


— CCq 1 

Ci(x) = — ^ — f r( y*) u g( x ,y*) (18) 

2U J 


where now 


r( y*) 


c x (y*) c(y*) 


-, 200 ^= 3 ^ 


U 



( 17 ) 


The only difference in evaluating y(y*) for horizontal and vertical 
gusts lies in the definition of the parameter having a variation of y* 
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across the span; for the vertical gust, that parameter is taken as the 
addi tional angle of attack and, for the horizontal gust, that parameter 

is 2oq — . The condition that 



7(y*) y* dy* = 2 


remains unchanged. 


Antisymmetric Span Influence Function y(y*) 

The antisymmetric span influence function 7 (y*) is defined over 
the span so that any given distribution of 7 (y*) will produce a unit 
rolling moment. These distributions refer to the span loading due to 
a linear angle of attack a = y* for the vertical gust or a linear 

leading-edge velocity 2oo — = y* for the horizontal gust. Four basic 

variations of y(y*) have been considered with the proper constants so 
that equation ( 15 ) is satisfied. The equations for the 7 (y*) varia- 
tions considered are given in table I and plots of these variations are 
shown in figure 3(a) . The names given to the four distributions obtained 
by rolling the wing refer to the distributions which would be produced 
by a uniform angle of attack. 


Correlation Function of the Rolling Moment 

Vertical gusts .- The autocorrelation function of the roll i ng moments 
due to vertical gusts at any two stations along the path of the wing is 
defined as 


k /.x Cl N ° l ( Xl ) 


(18) 
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With the substitution of the expressions for 0^x2) and C^(x;l) as 
given by equation (13), the correlation function of equation (18) becomes 


Sc ; ( x 2- x i) 



^ j i r(y 1*) w g ( x i^i*) d yi*/ l 7 (y 2 *) ^(*2^2*) ^2* 
r(y^)r(y 2 *) ^ / x v g( x i^i # ) » g (*a^ 2 *) ^1* ^2* 

r(yi*) r(y 2 *) ^2* 


(19) 


where it is assumed that the functions are convergent under either order 
of integration. An expression is thus obtained for the correlation func- 
tion of rolling moment in terms of the correlation function o±' vertical- 
gust velocity. In equation (19), 


$Wg( x 2-xi,y 2 *-yi*) = ^lim i J v g( x 2^2*) **1 

is the same as the two-dimensional correlation function defined earlier 
as equation (5) with x 2 - x-^ = Ac and y 2 * - yp* = Ay* = ^7— • 


By the proper substitution of variables, the double integrals of 
equation (19) may be separated into the single integral of the product 
of the integrated weighting functions of /(y*) and the correlation 

function i£w ■ Thus, with the substitution of 
8 


x 2 - x^ = Ut 



/ 
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equation ( 19 ) becomes 


where 



c i 2 r 2 

— 2- / r(n) t w (ut,ti) dTi 
8U 2 J 0 W g 


r(n) = 



^i*) ^i* 


and 


( 20 ) 


( 21 ) 


S Wg (uT,n) 



( 22 ) 


Equation (22) may be recognized as being equivalent to equation (5) • 

Horizontal gusts.- In an identical manner, the autocorrelation func- 
tion of rolling coefficient due to horizontal gusts at any two stations 
along the path of the wing is derived by use of equations (l6) and (l8) : 




ol rl 


/ 1 / 1 7 ( y i*) 7 ( y 2*) ^i* ^ 2 * 

With the same change of variables as in the preceding section. 


cto^Ci 2 p 2 

fci (Dr) - X W Ur ’ ,l) 


2U 2 


dn 


(23) 
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where 



(24) 


and equation (24) is now the equivalent of equation (4). The integral 
weighting function r ( tj ) is the same for both the horizontal- and the 
vertical-gust contributions to their rolling -moment correlation functions. 


Integral Weighting Function r(q) 

The integral weighting function r(q) as defined by equation (21) 
has been evaluated for the four distributions of y*) given in table I. 
These values are listed in table II and plotted against q in figure 3(b) . 
It may be shown that the nature of the function is such that the 
relationship 



r(q) dq = 0 


(25) 


must be satisfied for any variation of X .which pertains to an anti- 
symmetric variation of y(y*) . In table II the elliptic distribution 
is given in terms of K(k) and E(k), which are complete elliptic inte- 
grals of the first and second kind, respectively, of modulus k = ^ — — ■ 

The derivation of the elliptic weighting function is included in the 
appendix of the paper. 


Power Spectra of the Rolling Moment 

The power spectrum of the rolling-moment coefficient is defined 

as the Fourier transform of the autocorrelation function of : 
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® c 2 (^) = ^ J $c^(Ut) e u d(UT) 


( 26 ) 


For the vertical gusts , the power spectrum of the ro lling -moment 
coefficient may be found by substituting the derived relationship for 
(Ut) given by equation (20) into equation (26): 


$ Cz (cu) 




$ w _(Ut,ti) dTi d(Ut) 

O 



Changing the order of integration here is permissible inasmuch as the 
integrals of the correlation function of Vg are convergent in both Ut 

and q. The integral ^w is defined as 



P 00 _i“Ut 

/ e U $w I ,( UT > T l) 

J _co 6 


( 28 ) 


Similarly, the power spectrum of rolling-moment coefficient due to 
the horizontal component of gust is obtained from the substitution of 
equation ( 23 ) for the term fg (Ut) appearing in equation ( 26 ) : 

X 






an 


(29) 
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where the integral Iu in equation ( 29 ) is defined as 




-i— U t 


e U 


^Ug(UT,i]) d(Ut) 


(50) 


Thus, for two of the three components of turbulent gust velocities, 
the power spectrum of the rolling-moment coefficient is dependent on the 
integration of a function of the lifting distribution of the wing times 
a function which represents the Fourier transform of the correlation 
function of the vertical and horizontal gust components over the wing 
span. 


As previously stated, these results are based on quasi-steady con- 
siderations. Unsteady- lift effects can be taken into account simply by 
multiplying the power spectral density of the rolling moment due to each 


gust component by the function 



cp is the Sears function 


given in reference 10. 


Approximation for Side Gusts 

As pointed out previously, the side gust is treated here only in an 
approximate manner; that is, the spanwise effect is neglected. Based on 
this approximation, the rolling -moment coefficient is defined as 


C z (x) = c z 


Vg(Ax) 


U 

The correlation function is defined by 


(31) 


%v 


V UT) - — 

and the power spectrum is defined by 


s(Ur) 


(32) 


®c 


u 2 


G(o>) 


(33) 
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where 


gM 


_i_ 

3tU 



g(U-r) d(Ur) 


(3*0 


Relations Between the Yawing and the Rolling Moments 

No attempt is made herein to calculate directly the yawing moment 
due to atmospheric turbulence. Because of the more complicated nature 
of the phenomena which give rise to drag, as compared with those which 
give rise to lift, such an undertaking would be quite difficult. Further- 
more, in view of the fact that the yawing moments on the wing due to tur- 
bulence are relatively small, a detailed analysis would not generally 
be warranted. In this section, therefore, an approximate procedure is 
outlined for obtaining the yawing moments from the rolling moments. 

The yawing-moment coefficient due to sideslip can be expressed in 
the form 


On — Cnp(ct) p 


where. In this case, a is the sum of the trim angle 

instantaneous mean vertical-gust angle and where 

taneous mean side-gust angle so that 

TT * 


Oq and the 
p is the instan- 


Gp. = Cn^Ccfc,) 


Zs + (!%) Zs Zs 

u \ do, /„ U U 


where the second term is of higher order and is neglected. Similarly, 
differences in v g along the span give rise to higher order terms. 


The rolling moment can be expressed in the same form, so that the 
relationship between the yawing and roll trig moments due to side gusts 
is given by 


Gn( v g) 



c l( v g) 


(35) 
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Actually, this contribution to the yawing moment is generally negligible 
and is included here primarily for the sake of completeness. 

For the yawing moments due to vertical and horizontal gusts, similar 
reasoning may be employed. The yawing moment in these cases arises from 
the antisymmetric part of the instantaneous angle-of-attack distribution 
due to turbulence, as does the rolling moment, so that the two moments 
may be expected to be approximately proportional to each other; that is, 


^n( w g) 



( 56 ) 


Cn( u g) 



( 37 ) 


In essence, these relations imply that the yawing moment due to a given 
instantaneous spanwise gust distribution is jfche same as the yawing moment 
due to a linear gust distribution which has the same rolling moment. The 
deviation of the actual distribution from a linear one results in small 
differences in the vortex field and, thus, in small differences in the 
induced downwash. These differences lead to a contribution to the yawing 
moment which is believed to be small and, hence, has been ignored. 

In terms of their power spectra, the yawing moments are ^ defined as 


4, c n (“) 




"g 




Ug 


$ Cn(co) 






n P 


v g \Clf, 


'n T 


w g \° l P>cfc 


®C 7 (cd) 


'6 




w, 


g 


(38) 


The power spectra of the rolling moments are defined in the preceding 
sections . 
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APPLICATION 


Approximations to the One-Dimensional (Point) 

Correlation Functions 

In order to evaluate the effects considered in the preceding part 
of the paper, calculations will now be made by using the results derived 
therein. These calculations will be based on a simple analytical expres- 
sion for the longitudinal point correlation function which has been sug- 

of measurements in wind tunnels: 

_ki 

e L (39) 

of turbulence defined for any longl- 

by 

f(r) dr (40) 


The characteristics of clear-air turbulence measured in the atmosphere 
(ref. 13 ) may be shown to be reasonably well represented by equation (39 ) , 
with a value of L of approximately 1,000 to 2,000 feet. There are some 
theoretical objections to this function - primarily the fact that it has 
a nonvanishing slope an r — =>0 and, hence, that the associated power 
spectrum does not decrease rapidly enough for very short wavelengths. 

These conditions imply that the mean square of the derivative of the gust 
velocity with respect to the space coordinate is infinite. However, from 
available measurements on atmospheric turbulence, it appears that equa- 
tion (39) remains valid to distances which are small compared with the 
span of the airplane (on the order of several inches), and the behavior 
of the spectrum at very short wavelengths is relatively unimportant because 
airplanes c ann ot respond to them to any appreciable extent. Therefore, in 
the absence of more reliable information all calculations described in 
this paper are based on equation (39)- 

The corresponding lateral correlation function related to f(r) by 
equation (2) is found to be 


gested in reference 12 on the basis 

f(r) = 

where L is the longitudinal scale 
tudinal correlation function f(r) 
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A plot of the functions given by equations (39) and- (4l) is shown 
in figure 4. Their respective power spectra, denoted by G(k') and 

F(k') where k' = — , are- given by 


F(k') = S 

rtU 


1 + (k')‘ 


( 42 ) 


G(k‘) = ±- 
jtU 


1 + 5(k') 2 
[l + (k') 2 ] 2 


(43) 



Calculations for Vertical Gusts 

Rolling-moment correlation function. - When eq uation (4l) is subBti- 
tuted into equation (22) with r = y( Ur) 2 + the correlation func- 
tion defined by equation (22) becomes 


$ Wg (UT,T]) = w g 2 


1 - 



(44) 


Inasmuch as the evaluation of the rolling-moment correlation function, 
as such, is not necessary to the analysis of this paper, only limited 
consideration is given to the calculation of autocorrelation functions. 
Equation (20) has been evaluated in closed form for the case of the 
rectangular distribution of the span influence function y(y*) as given 
in tables I and II: 
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where Kq and Kq are defined in reference 6 an incomplete modified 
Bessel functions where 


and 


, ,_iP 

^ o sinh x — 

K v (p',A) - / * e'* cosil 6 cosh v8 d0 (46) 

J 0 



These two parameters represent the ratios of the distances b and Ur 
to the integral scale of turbulence L. The parameter p 1 reflects the 
size of the wing span relative to the characteristic size of the turbu- 
lence and, as such, is one of the more important parameters appearing in 
all the calculations involving spanwise correlation. It effectively 
scales the magnitude and shape of the correlation functions and power 
spectra and, in the limit as p ’ — >0, the equations for the antisymmetric 
moments likewise go to zero inasmuch as no rolling or yawing moment will 
exist when a finite span shrinks to a point. 

The parameter A is a measure of the flight-path distance relative 
to the characteristic size of the turbulence and, in the limit as A — >0, 
the correlation function must reduce to the mean square value of the 
rolling-moment coefficient; hence, 

G? - $Ci(* = 0) 


>g 2 Cz p 2 
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+ 12p’ 2 


+ 24p' + 24) e”^ 
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(47) 



NACA TN 3864 


24 


Inasmuch as no adequate tables appear to be available for the func- 

rv ~ . 

tions Kq and Kq, a numerical evaluation of equation (45) has not been 

made. Hove ver , an analysis of this correlation f unction vith other 
approximations for f(r) indicates that the effect of span loading is 
minor and that a reduction in (3 1 attenuates the correlation function". 


Evaluation of 



- For the vertical gust component. 


the inte- 


gral definition of is given by equation (28) whereas (Ut,t}) 

is now defined by equation (44). The indicated integration may be per- 
formed in closed form as a function of rj and the reduced frequency 
parameter k' . Thus, 



where 



and Kq and Kq are modified Bessel functions of the second kind of 
order 0 and 1, respectively. 

A plot of equation (48) is shown in figure 6 as a function of the 
frequency parameter k', for a range of values of p'q/2 from 0 to 1.0. 
Although the physical significance of the function I w is rather obscure, 
the plots are useful in the subsequent numerical integration of the prod- 
uct of Iw and F. 

Power spectrum of rolling moment .- In general, the analytical solu- 
tion of equation (27) for the power spectrum of the rolling-moment coef- 
ficient due to vertical gusts, when possible, is a tedious process. 
Numerical integration by means of either SimpBon’s rule or some inte- 
gration process of higher order is generally preferable to integration 
in closed form. However, the analytical evaluation of equation (27) for 
the case of a wing with rectangular span loading is given here in order 
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"bo illustrate some of the characteristics of the equations. After the 
indicated substitutions are made, equation ( 27 ) becomes 




8jtu3 


J (4 - 6 t! + n 5 ) 




P 


2 / 2 

7 Kq + 


^r 1+ 

? L 


5(k')‘ 


1 + ( k ') 2 £ + ( k *) 2 ] 


5/2 


Kn 


dq 


for which the integrated solution is 




- — -(a3(k*) 2 P Ko(x) cbc + ia^ + l6a 2 ("l - (k') 2 ]l K<j(a) + 

(k-) 2 ]\ ° l L J J 

Isb?\3 - (k') 2 ] + 32a [l - (k’) 2 ]j K x (a) + 2a 2 [l - 3(k’) 2 ] - 32 [l - (k') 2 ] j (4-9) 


where a = p 'Ul + (k') 2 , k’ = — , and Ko(a) and K-i(a) are modified 

* U 

Bessel functions of the second kind of argument a.^ Equation (49) is 
plotted in figure 7 (a) as a function of k’ for a range of p 1 between 
0.03125 and 1 . 0 . 

For small values of frequency a> (and hence k 1 ) or scale factor p ' , 
equation (49) becomes poorly behaved because the solution takes the form 
of small differences of high-order terms. The reason for this may be seen 
by expanding the Bessel functions in their power-series form a.nri grouping 
like powers of the variable a. The coefficients of the first three terms 

of the power series a a - ^, and (which are the predominant terms 

for values of a < 1) are identically zero. Older these conditions, small 
computing errors or the lack of significant figures will cause large inac- 
curacies in the numerical evaluation of the function. 


The difficulties just described may be overcome somewhat by evaluating 
equation (49) for the limiting case of k' = 0 : 


® c ^(k'= 0 ) = 


l 8 w g 2 LC z 


nuV* 


(p' 1 * - + 16 p’ 2 )ko(p') + ^6p’ 5 -s- 32p *) KjlCp 1 ) + 2P’ 2 - 32 

( 50 ) 


^Values for the integral of Kq may be found in several publica- 
tions, one of which is reference 14, table 2 (Zahlentafel 2). A com- 
prehensive listing of other available mathematical tables including these 
Bessel functions is given in reference 15. 
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When the Bessel functions are again expanded in powers of p * , only several 
terms are needed to evaluate the function at small values of p'. As before, 
the coefficients of all negative orders and the zero order of p 1 are 
identically zero. 

The physical necessity that, as the span b approaches zero, the 
expression for the power spectrum of the rolling-moment coefficient must 
also approach zero is satisfied by equation (49) inasmuch as the lowest 
order term with a nonzero coefficient appe_aring in the equation is a 2 
(as pointed out above); that is, for b — >0, 


®Q^(a>) ** (Constant) a 2 = 0 


In order to compute for the other three types of distribution 

L 

of wing loading given in tables I and II, a numerical-integration process 
involving Simpson's three-point rule of integration was employed. The 
power spectra thus obtained are plotted in figures 7(b), (c), and (d). 

This method was also used for the rectangular lift distribution and was 
found to give good agreement with the analytical results. 

It is of interest to note that whereas the power Bpectra of the 
vertical gust approach a logarithmic decrement of -2 (see fig. 5 ) > the 
rolling-moment power spectra shown in figure J approach a decrement of -J. 
At the low-frequency end of the spectrum (long wavelengths) the power 
appears to approach a constant which is zero only when p ' , the ratio of 
span to scale of turbulence, is zero. 

Some simplified approaches to the calculation of the rolling power 
of gusts (for example, ref. 4) lead to the result that the spectrum of 
the rolling power of the vertical gust appears as the first derivative 
(slope) of the vertical-gust spectrum. As may be seen from figure "J, 
such an approximation is justified only in a very small band of frequencies 
for wings having small values of p ' . 


Calculations for Horizontal Gusts 

Rolling-moment correlation function .- When the expressions for f(r) 
and g(r) given by equations (39) and (4l) are substituted into equa- 
tion (24) with r = 
tion for horizontal gusts becomes 


I 


(Ur ) 2 + 


bry 

2 


the one -dimensional correlation func- 
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The correlation function of rolling moment is obtained by inserting equa- 
tion (51) into equation ( 23 ) and integrating. For a rectangular distri- 
bution of ^(y*). 
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(52) 


and 



48u g 2 a 0 2 C2 


uV* 


5 + 3P' 2 + 6p' + 




(53) 


As in the case Of the vertical gust, the correlation function has not 
been calculated for the other three distributions of /( y*) for the 
reasons alr eady given. 

Evaluation of I u ^, ^.- The evaluation of I u , as defined by equa- 
tion ( 30 ) for the case of the horizontal gust, is given by the expression 
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where k' = — , and. Kq and 
the second kind of argument 


are the modified Bessel f unctions of 



The function given by equation (5k) is plotted against k 1 in fig- 
ure 8 for values of p'n/2. Values taken from this plot may be used in 
the numerical integration of the power spectra of the rolling-moment 
coefficient . 


Power spectrum of rolling moment . - The power spectrum of the rolling 
moment due to the horizontal components of turbulence acting on the wing 
has been determined by using the expression for I u obtained in the 
preceding section and the four distributions of the parameter r( ri) given 
in table IX. The integral of equation (29) has been evaluated numerically 
for all four cases of load distribution, and the resulting variations of 
the power spectrum with frequency and g ' are plotted in figure 9- In- 
addition, the analytical solutions for the cases of rectangular and par- 
abolic distributions are given here and their numerical values were 
checked against those obtained by the numerical- integration process. By 
use of equation (5k), ike solution for the rectangular case is found to 
be 
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l 6 a(a 2 + 4)K;j_(a) - 64 

where a = p* ]jl + (k 1 ) 2 and k 1 = The analytical solution for the 

parabolic distribution is given by 



240ug2La o 2 C2 2 r 

- (209a6 + 25,344a4 + 331, 776a 2 ) Kq (a) + 

7itU 3 a8 1 + (k 1 ) 2 ]*- ' 

(aT + 3,424a5 + 133 , 632 a? + 663,552a)Ki(a) + 

(aT + 63a?) J Kq(x) dx - 1,120a 14 " + 32,236a 2 - 663,552 


(56) 

A comparison of the values obtained for equations ( 55 ) and ( 56 ) and 
plots of the results obtained by the numerical- integration process indi- 
cated no difference, and none is shown in figure 9 . 

It is significant to observe that very little variation exists in 
the power spectra of figure 9 tor the four span loadings consi der ed. 
However, as compared with the rolling moment due to vertical gusts 
(fig. 7)> the rolling moment due to horizontal gusts is relatively small 
for small values of trim angle of attack. Although no exact expression 
for the ratio of the power spectra of the rolling moments due to Ug 

and Wg may be given without including { 3 1 and y, it may be seen from 
figures 7 and 9 that, in general, 
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u s 




0.2 



(57) 


where 


cxq is given in radians. 


Calculations for Side Gusts 

For the side gust considered, the correlation function of the 
rolling-moment coefficient as given by equation ( 32 ) becomes 



(58) 


and the variation of this function with Ut/l is, of course, equal to 
the variation of equation (4l) with r, which is plotted in figure 4. 


The power spectrum of the rolling -moment coefficient as given by 
equation (33) with G(k‘) given by equation (43) becomes 


%(*’) 


Ip V S L 1 + 5(k')‘ 


itu5 


_ / 

l + (k') 2 4 
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(59) 


The variation of the spectrum with frequency k 1 = is shown as the 
G(k') curve of figure 5* 


DISCISSION 


The purpose of this section is to discuss the implications of the 
assumptions made in the analysis of this paper, the reasons for making 
these assumptions, and the application of the results. 
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Assumptions Concerning the Nature of Turbulence 

The turbulence was assumed to be homogeneous in order to make the 
problem stationary in the statistical sense and thus permit the use of 
the mathematical techniques developed for such problems. In a practical 
sense, turbulence can be homogeneous only in a limited body of air. The 
assumption thus implies that the dimension of this body of air along the 
flight path is large compared with the distance traversed in the reaction 
time of the airplane. In the case of loads studies this reaction time 
is of the order of the time to damp to one -half amplitude, but, in the 
case of motion studies, the reaction time may be much larger. Obviously, 
the greater the body of air, the greater the reliability with which the 
loads and motions can be predicted (in a statistical sense) for one run 
through it. In general, turbulence at very low altitudes, which may be 
influenced significantly by the configuration of the ground, and the 
turbulence in thunderstorms may not be sufficiently homogeneous for this 
type of analysis, but other types of turbulence are likely to be sub- 
stantially homogeneous over sufficiently large distances. 

Isotropy was assumed in order to permit the required two-dimensional 
correlation functions to be expressed simply in terms of the one- 
dimensional correlation functions. For sufficiently short wavelengths 
all turbulence is isotropic (see ref. 7 ) > but for long wavelengths it 
can be isotropic only if it is homogeneous (both in the plane of the 
flight path and perpendicular to it). (The condition of axisymmetry is 
less restrictive inasmuch as it does not specify the variation of the 
characteristics of the turbulence in the vertical direction. ) In prac- 
tical problems, if the turbulence may be assumed to be homogeneous, the 
conditions of isotropy are likely to be satisfied sufficiently to permit 
the use of the approach presented herein for all but very long wavelengths. 
The wavelength at which this approach ceases to be valid depends on the 
size of the body of air under consideration, being longer for a large 
body. 


Taylor's hypothesis implies that the variation in gust intensity 
that prevails along the flight path at any instant will remain substan- 
tially the same until the airplane has traversed the given body of air. 
The required correlation functions for atmospheric turbulence are thus 
in the nature of space correlation functions (rather than time correla- 
tion functions) and have been considered as such. The statistical char- 
acteristics of the turbulence are then independent of the speed at which 
it is -traversed. Clearly, the validity of this hypothesis depends on 
the flying speed of the airplane and it would be expected that, at very 
low speeds, the hypothesis of Taylor becomes less valid and the results 
may be less accurate. On the basis of present knowledge, no definite 
lower limiting speed can be quoted. The effect of finite flying speed 
on the gust correlation function can be expected to be most pronounced 
for large distances, where the correlation is weak. Thus, the effect 
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on the various spectra is likely to be small and to occur at the longest 
wavelengths, where, as previously mentioned, the spectra are somewhat 
■uncertain for .other reasons as well. 

For practical purposes, the parameter L (the integral scale of 
turbulence) used herein is a largely fictitious quantity, inasmuch as 
it is, to a large extent, proportional to the -values of the gust spectra 
for infinite wavelengths . In view of the uncertainties in the values 
of the spectra at long wavelengths and the fact that the spectra in this 
region predominantly define the area under the integral, the parameter L 
has little physical significance. Therefore, at present, insufficient 
information is available to give an exact value for L to be used in 
connection with the numerical results calculated herein. However, on 
the basis of the measurements such as those of reference 13 , a value 
of 1,000 to 2,000 feet appears to be appropriate for the conditions of 
the referenced tests. It is desirable to obtain more inf ormation con- 
cerning the spectra of atmospheric turbulence under a wider range of 
conditions. More definite values could then be deduced by fitting meas- 
ured results by means of an analytical expression of the type used here. 
This expression could be used as a means of obtaining a value of L by 
extrapolation of the measured results to infinite wavelengths (zero 
frequency) . 


Assumptions Concerning the Aerodynamic Forces 

The fundamental assumption concerning the aerodynamic forces is 
that they vary linearly with gust intensity. This assumption implies 
that the ratio of the gust speed to the flying speed must always be 
fairly small; if the aerodynamic forces and moments tend to vary with 
gust intensity in a nonlinear manner, as the wing yawing moments do for 
all angles of attack and the other forces and moments do for high angles 
of attack, the ratio of gust intensity to flying speed must be very 
small - about 1/30 or less. However, as previously mentioned, the wing 
yawing moments due to gusts are likely to be quite small, so that some 
error in them due to slight deviations from linearity is not likely to 
affect appreciably the results of an analysis of the lateral motion. 

Hence, for an airplane flying at small angles of attack and at speeds 
of about 200 knots or more, in continuous turbulence, the assumption of 
linearity should be valid; for flight in severe thunderstorms, it Is 
not likely to be valid, and, for flight at high angles of attack, it is 
likely to be valid only for light turbulence. 

The rigidity of the wing, which was mentioned in the list of assump- 
tions, enters only indirectly into the problem considered herein. The 
results obtained here are valid whether the wing is rigid or not. How- 
ever, in the case of flexible wings (the term "flexible" being used to 
describe-'' wings with deformations which give rise to appreciable aero- 
dynamic forces), certain additional information is required. (See ref. 6 .) 
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This information may take the form of span influence functions y(y) 
modified by static aeroelastic effects, or may require certain cross- 
correlation functions or cross spectra between the gust forces and the 
dynamic forces, depending on the individual case. 

The assumption that the indicial-response influence function h(t,y) 
can be written as a product of functions of time only and distance along 
the span only is based on the reasoning of reference 6 . This reasoning, 
in turn, is based on the observation that, according to the available 
information for the lift distributions due to sinusoidal motions (and, 
hence, those due to indi cia] motion), the lift distribution tends to be 
substantially invariant with frequency (or time) except for an overall 
factor. Inasmuch as this information is confined to unswept wings, this 
assumption may not be valid for swept wings. 


Application of the Results 

In this paper the rolling moments and yawing moments have been cal- 
culated for a wing due to the u, v, and w components of turbulence. 

If the turbulence is isotropic, these components are statistically inde- 
pendent at a point. In any practical application, all three components 
are always present and the wing rolling and yawing moments due to the 

combined action of the three components must be known. In isotropic 

turbulence, the cross correlations between u and w and between v 
and w in the horizontal plane are zero, although u and v have a 
nonvanishing cross correlation. Thus, the moments due to v and w can 

be added directly, but, if horizontal-gust effects are to be taken into 

account, not only the moments due to u calculated herein but also the 
moments which arise from the cross correlation between u and v should 
be added to the others. However, there is reason to believe that the 
horizontal-gust effects on the lateral moments are generally very small, 
so that neglect of this cross -correlation effect is usually justified. 

The rolling and yawing moments due to Ug and Wg considered herein 

are only those contributed by the wing but, inasmuch as the lateral moments 
contributed by the fuselage and tail as a consequence of these two com- 
ponents of gusts are generally very small, the results given here may, 
in general, be used to represent the lateral moments on a complete air- 
plane due to these two gust components. 

Similarly, the rolling and yawing moments of a complete airplane due 
to the v-component of gusts depend not only on the wing contribution 
considered here but also on the contribution of the vertical tail, 
which can be calculated in a straightforward manner. For instance, a 
method of calculating the yawing moments and side force on a fuselage 
and vertical fin due to side gusts is found in reference 17 . 
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Although the contribution of the horizontal component of gusts to 
the lateral moments appears to be small compared with the other two com- 
ponents, it should be kept in mind that the effect of this component 
increases as the square of the trim angle of attack. (See eq. (57)0 
For conventional airplanes in the landing configuration and for verti- 
cally rising airplanes in the transitional stage, the effects of hori- 
zontal gusts may well be predominant in calculations of the forceB, 
moments, and motions due to turbulence . 


CONCLUDING REMARKS 


The correlation functions and power spectra of the rolling and yawing 
moments on an airplane wing due to the three components of continuous 
random turbulence have been calculated. The rolling moments due to the 
longitudinal (horizontal) and normal (vertical) components depend on the 
spanwise distributions of instantaneous gust intensity, which were taken 
into account by using the inherent properties of symmetry of isotropic 
turbulence . The results consist of expressions for the correlation func- 
tions and spectra of the rolling moment in terms of the point correlation 
f -unctions of the two components of turbulence. 

Specific numerical calculations were made for a pair of correlation 
functions given by simple analytic expressions, which fit available 
experimental data very well. Calculations were made for four lift dis- 
tributions and the differences in the results calculated for these dis- 
tributions were small. By comparison with the results calculated herein, 
the results of previous analyses for which it was assumed that random 
turbulence along the flight path and variations of turbulence across the 
span were linear have been shown to be valid only when the ratio of the 
span to the integral scale of turbulence (about 1,000 to 2,000 feet) is 
small. 


A comparison of the power spectra of the rolling moments due to 
horizontal gusts and those due to vertical gusts showed that the vertical 
gusts were predominant at small values of trim angle of attack (or trim 
lift coefficient) ; however, the relative effect due to horizontal gusts 
increased as a function of the square of the trim angle of attack. 

The rolling moment due to lateral (side) gusts, which is small, was 
expressed in terms of the instantaneous value of the gust at representa- 
tive points on the wing, so that the effect of spanwise variation in gust 
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intensity was ignored. The yawing moments were considered to be propor- 
tional to the rolling moments, the constants of proportionality being 
given by simple aerodynamic relations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 6, 1956. 
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APPENDIX 


EVALUATION OF THE ELLIPTIC INTEGRAL WEIGHTING FUNCTION 


The evaluation of the integral weighting function r(Tj) involves 
the integral given by equation (21) : 


F(t]) = J r(yi*) 7(y!*+Tl) dy-j* 

For the case of the elliptic distribution of the additional span loading 
factor. 


r(y*) = ^ y*\Ji - y * 2 

and the integral, weighting function to be evaluated becomes 
2 1 TJ 

r( ’ 1) ‘ (¥) /_ x (yi*)(yi* + 1 *) V 1 - yi* 2 ^ 1 - ( y i* + 1 ') 2 

Older the substitution 


2y 1 * + H 2y 1 * + T] 

x = ■ — 

2 - TJ 6 


the integral may be written as 

r(n) = (if) (fir) f_ ± (8 * " + - ti) 2 ^1 - i(5x + nf dx 
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Inasmuch as 1 - z 2 = (l - z) (l + z ) , 


rtn) - (f ) (H 21 ) ( 82x2 - 1 ' 2 )^/( 1 - 1 8x + 2 ”) ( x + 2 8:5 " 2 
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With the notation k = = — - — L, 

2 + T] 2 + 1} 

r(n) = ^|(2 - ti) 2 (2 + Ti) J |(2 - T }) 2 * 2 - i} 2 ] \ji 1 - kx) (l + x) (l - x)(l + kx) dx 

3t — X 

= ^(2 - T!) 2 (2 + n) ^[(2 - n) 2 x 2 - Ti 2 ]^(l - k2x2)(l - x 2 ) dx 


where the integrand may be seen to be an even function of the variable x. 
Multiplying numerator and denominator by the radical and expanding yields 


r(n) - %2 - n) 2 ^ + n) 

XT' 


*2(2 - n) 2 r 1 T - »®[ca + n) 2 + (2 - *0 2 + n 2 ] f 1 , ^ 

l/(l - k2»2)(l - *2) L JJ 0 : k2x 2)(! . “ 


[( 2 - n) 2 + n 2 + n^] J 1 ■. ** f 1 25 

0 p. - k2x2) (l . x 2) J 0 |j(i . k 2 x 2)(i - *2) 


The integrals may be recognized as elliptic integrals in powers of x 211 
for which the closed-form -solutions may be found in reference 16 , for 
example. In terms of the standard elliptic integrals (in Jacobi's nota- 
tion) of modulus k = 2 — H . 

2 + n 


dx 


1 _ x 2) (1 - k2 x 2) 
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•which is defined, as a complete elliptic integral of the first kind, and 

. . r 1 / 1 - k 2 x 2 

E(k) = / l/ dx 

J 0 f 1 - x 2 

•which is defined as a complete elliptic integral of the second kind. 
Tables of these integrals may be found in most mathematical handbooks 
as well as in reference 1 6 . In terms of these integrals, the solution 
for the integral weighting function is found to be 


r(t|) => ^<2 - T!) 2 (2 + n) \ 

rt 2 


k2(2 - tj)‘ 


15 k P 


(8 + 3k 2 + Ik^)K(k) - (8 + 7k 2 + Sk^ECk) 


k 2[(2 + ,) 2 + (2 - n ) 2 + d n + k2 ) K(k) . 2(l + k 2 ) E(k) 


3k 


[(2 - n ) 2 + n a „* Jgj[- K(k) . E(k) j . ,% (k) 
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TABLE I 

VARIATION OF y(y*) 


Distribution 

r(y*) 

Rectangular 

6y* 


Elliptic 


y*2 

Parabolic 

i5y*(i - 

y* 2 ) 

Triangular 

24y*(l - 

M) 


TABLE II 

INTEGRAL WEIGHTING FUNCTION r(t]) 


Distribution 

r(n) 

Limits 

Rectangular 

6(4 - 6rj + t]3) 


Elliptic 

2 + t]) [4t](ti 2 _ 3 n - l)K(k) + 

15* 2 

(4 + 9n 2 - H^)E(k)j 

0 = T) — 2 

Parabolic 

||(64 - 336t[ 2 + 280t)3 - 42^5 + 5ij7) 

0 = i] = 2 


^^(2 - 10T| 2 + 5ti3 + 5^ - 3tj5) 

15 

mm 

Triangular 

m 

—(8 - 20ti + 10t} 2 + 5t]3 _ 5^4 + Tp) 

15 
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(a) Wing passing through three-dimensional turbulence. 


Longitudinal components, f(r) 




Lateral components , g (r) 


(b) Components of turbulence as a function of distance r. 

Figure 1.- Sign convention and stability axes of a wing passing through 
a turbulent velocity field. Arrows denote positive direction, where 
applicable . 
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(b) Side gust components. 

Figure 2.- Schematic drawing of the relationship between the components 
of horizontal and side gusts at any two arbitrary points. 




fa) Variation of span influence function. 


(b) Variation of integral weighting 
function. 


T34 *2 .4 ^4-4 ^-P 4 _ -P’1 4?,, 4-4 -w .3 4 ~ 4-^ 1 T «^ rt V4-^/v 

J? i^LLL C ^ • — f OJIJLO. ljA,\JLL U1 DJJCULL ±1LLJ,UC1UJC J. UJU.^ ^-LWIX / CLLLU. J-LLUC^J, CM. WgXftm/XUfe 

function T for four types of loading. 
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Figure 5-~ Power spectra of lateral and longitudinal components of 
isotropic atmospheric turbulence . 
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Figure 6.- Variation of vertical-gust weighting par ame ter 


Iw(*') 


for 


a number of values of 3 t]/2. 


kd 
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(a) Rectangular span loading. 

Figure 7*- Power spectra of rolling moment of wing due to vertical gusts. 
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(b) Elliptic span, loading . 


Figure 7*- Continued. 
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Figure J.- Continued. 
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Figure 8.- Variation of horizontal-gust weighting parameter 


IuU') 

u s 2l 


for 


a range of values of 0t]/2. 
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(b) Elliptic span loading. 

Figure 9* - Continued. 
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Figure 9 • - Continued. 
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NACA - Langley Field, Va. 


